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ABSTRACT 

A brief description of the library of the Standard Model and anomalous WWj coupling 
contribution to the matrix element for e'^e~ uunj process is given. It can be used 
with any Monte Carlo program for e~^e~ —>■ ffwy processes. A working example of the 
application for the KORALZ version 4.04 is also provided. 
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1 Introduction 



The z^z>7 production in e^e~ collisions is of great interest, as it is sensitive to the triple 
gauge boson coupling WW'y of the Standard Model. Its precise measurement will not 
serve only as a stringent test of the SM, but may reveal (or constrain) anomalous gauge 
couplings. The events with photon(s) plus missing energy might originate also from 
other mechanisms, signalling new physics beyond the Standard Model. For example, 
such final states can be produced in both gravity- and gauge-mediated supersymmetric 
models or low-scale gravity. The missing energy in these events is caused by weakly 
interacting supersymmetric particles, such as gravitinos, neutralinos and/or sneutrinos 
or gravitons 0]. In all such cases the Standard Model e~^e~ uu'y events are irreducible 
background and reliable theoretical predictions for them are therefore necessary. 

With the sensitivity afforded by the LEP experiments [0, |], and expected at future 
e'^e~ colliders 0, the photon(s) plus missing energy events provide an opportunity to 
search for new physics phenomena. Any meaningful interpretation of the experimental 
data requires a Monte Carlo simulation in which Standard Model predictions may be 
augmented by the contributions from possible anomalous couplings. 

Since the LEP collaborations are entering their final years of operation, now is a good 
time to document the programs that have actually been used in the data analyses at 
LEP. In the present paper we document the library for the calculation of the effects of 
the WW'y interaction for the e^e~ —>■ vvn'y process within the Standard Model as well 
as from the anomalous couplings. It is based on the work of ^ and the description of the 
physical content of the program interface and discussion of its uncertainties can be found 
in 1^. For an alternative implementation of the WW'y vertex, see 

In principle our library can be combined with any e+e" ff ri'y generator, but in 
the present paper we will use an interface to KORALZ version 4.04, described in []TU|, ITT 



as a working example. That is the reason why the fortran code of the library will be 
archived together with the KORALZ tree of subdirectories 0]. Let us note that, in 
future, KORALZ will be replaced by a new program, KK2f |T^, which is based on a more 
powerful exponentiation at the spin amplitude level and in which the library will also be 
easy to implement. 

Version 4.04 of the KORALZ Monte Carlo program can be used to simulate e^e~ — 
// n7, (/ = fi, T, u, d, c, s, b, u) processes up to the LEP2 energy range, including YFS 
exclusive exponentiation of initial- and final-state bremsstrahlung and, optionally, the 
effects of various anomalous couplings. In the case of the LEP2 centre-of-mass energies 
and if f = u, the present library can be used for that purpose. 



2 Calculation of anomalous couplings 

To evaluate the effects of anomalous WW'y, a tree-level calculation of the squared matrix 
element for the process e+e" lyuy has been carried out 0. It includes the effects 
of the anomalous C- and P-conservingQ contributions parametrized with the help of the 

^For a recent phenomenological analysis of CP-violating couplings in e+e^ vv^, see 0. 
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couplings Ak^, A-^ (in what follows we will suppress the subscript 7). The library is 
formed on this basis. When activated, it uses the 4-momenta of the neutrinos and the 
photon provided by the host program to compute a weight, w, for each event according 
to 

I jlaWW-y excl. . A^WWy . kAWW'y\2 
^ K^'SM +-^^'SM +-^^'ano I /-■ \ 

I A/f^^T ea'c/.|2 • 
K^'SM I 

M^^^ is the matrix element due to the anomalous 7^ 0, A 7^ couplings, the A^sm^''' 
is the matrix element due to the SM WWj interaction and A^sm^'*^ represents the 
remaining part of the SM matrix element for the e~^e^ uu'-f process^. Note that because 
of the above separation, even for the Standard Model WW'^ interaction, the use of the 
library is necessary to calculate the weight w. 

As the calculation of ref. is performed at 0{a), the case of multiple bremsstrahlung 
requires a special treatment. In this case, a reduction procedure is first applied in which 
all photons except the one with the highest-p^ are incorporated into the 4-momenta of 
effective beams. In the second step, the 4-momenta of the highest pt photon, the effective 
beams and neutrinos are then used to compute the weight. Cross-checks of the calculation 
as well as checks of the validity of the reduction procedure are described in . The results 
of the calculation have been used in the measurement of the WW'j coupling described 
in m. 



3 Flags to control anomalous couplings in KORALZ 

The calculation of weights for anomalous couplings is activated by setting the card 
IFKALIN=1. This is transmitted from the main program via the KORALZ input parame- 
ter NPAR(15). Additional input parameters are set in the routine kzphynew(XPAR,NPAR), 
although there are currently no connections to the KORALZ matrix input parameters 
XPAR and NPAR Q Table |I| summarizes the functions of these input parameters. 

More input parameters are initialized in the subroutine anomini, which is placed in 
the file gengface.f and subroutine initialize of the file geng.f. Both files are placed 
in the directory korz_new/nunulib. 

In order to provide the user with enough information to retrieve the w for a given 
event for any Ak, A, we take advantage of the fact that, for each event, one may write the 
It; as a quadratic function of the anomalous couplings, in terms of the results calculated 
for six numerically distinct combinations of the Ak, A values as follows: 

A 2 / x2 ^ An \ / Ak 1 Ak 2 



MAk,A) = (l-(-r-(_)^ + -_)^0,0)-(— --(— )^)„,(-Ak„,0) 

/ Ak I.Ak.o a Ak \ ,^ ^, /A 1/A,2\ ^ N 



^Note that such a separation is gauge-dependent and, if not treated carefully, could lead to meaningless 
results. See ref. for details. 

^ In most uses of KORALZ, the numerical value of these parameters is irrelevant or defaults are 
sufficient. It is expected that the advanced user may like to change them, connecting directly the 
kzphynew (XPAR, NPAR) routine with her or his main program. 
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Parameter 


Description 


Default 


lENRICH 


Enrich spectrum of generated sample with hard-photon events 
IENRICH= 1/0 (on/off) 





IRECSOFT 


Generate only events with photon(sj ii IREubUr 1=1 





EMINACT 


Minimum sum of all photon energies required to calculate 
anomalous weights 


17 GeV 


EMAXACT 


Maximum sum of all photon energies allowed to calculate 
anomalous weights 


1000 GeV 


PTACT 


Minimum sum of all photon momenta transverse to the beam 
direction required to calculate anomalous weights 


2 GeV 



Table 1: Input parameters to control the calculation of weights for anomalous WW'y 
couplings. 



/A 1,A,2 A Ak \ , , , A Ak , s , . 

When the calculation is completed, the six weights are stored in the common block common 
/kalinout/ wtkal(6), with the assignments shown in Table The user is then free to 



Common block entry 


Weight parameter 


wtkal(l) 


w{0,0) 


wtkal(2) 


w{—Ako, 0) 


wtakl(3) 


w{Akq, 0) 


wtkal (4) 


wiO, -Ao) 


wtkal(5) 


w{0,Xo) 


wtkal (6) 


w{Ako,Xo) 



Table 2: Correspondence between entries in kalinout common block entries and weight 
parameters of eq. |^. 

calculate the w for whatever combination of Ak and A is desired. In our program we set 
Akq = 10 and Ao = 10. If the lENRICH input parameter is set to 1 the generated sample 
will have more events with hard photons than predicted by the Standard Model. The 
appropriate compensating factor is included into the weights wtkal. It is thus always 
assured that e.g. the generated sample, if the weight wtkal (1) is used, represents the 
Standard Model predictions. 

The code for a calculation of the weight w is placed in the directory korz_new/nunulib 
in the files geng.f and gengface.f. 
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4 Demonstration programs 



The demonstration program DEMOS . f for the run of KORALZ when our hbrary is acti- 
vated can be found in the directory korz_new/f ebruary and the output DEMOS. out in 
the directory korz_new/f ebruary/prodl. The DEMO.f for the run of KORALZ with the 
Standard Model interactions only and its output DEMO . out are also included in the direc- 
tories mentioned above. All these files as well as the library itself are archived together 
with KORALZ [O]. 
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